Ultracold three-body collisions near narrow Feshbach resonances 
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We study ultracold three-body collisions of bosons and fermions when the interatomic interac- 
tion is tuned near a narrow Feshbach resonance. We show that the width of the resonance has a 
substantial impact on the collisional properties of ultracold gases in the strongly interacting regime. 
We obtain numerical and analytical results that allow us to identify universal features related to 
the resonance width. For narrow resonances, we have found a suppression of all inelastic processes 
in boson systems leading to deeply bound states and an enhancement for fermion systems. 
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Strongly interacting three-body systems play an im- 
portant role in a great many areas of physics includ- 
ing condensed matter, atomic, molecular and nuclear 
physics [1, 2]. With the recent advances in the control 
of interatomic interactions, ultracold atomic gases have 
rapidly become a preferred test bed for many interesting 
physical phenomena. One of the most important tools 
for this control is the Feshbach resonance [3]. Applying 
a magnetic field, the s-wave two-body scattering length 
a between two atoms can be tuned from — oo to +oo, al- 
lowing for the study of systems from non-interacting to 
strongly interacting. A major experimental difficulty en- 
countered in the strongly interacting limit, \a\:^ro where 
tq is the characteristic range of the interatomic interac- 
tions, is that the system becomes unstable due to three- 
body collisional loss of atoms and molecules. 

Ultracold three-body collisions in the limit |a|;»ro 
are intimately connected to a fundamental and counter- 
intuitive effect first predicted by the nuclear physicist Vi- 
taly Efimov in the early 1970's [4]. He showed that three- 
boson systems with |a|^oo possess an infinite number of 
weakly bound three-body states even when there are no 
two-body bound states. These Efimov states have a pro- 
found effect on low-energy scattering observables as they 
lead to interference and resonance effects [6, 7] that pro- 
duce distinctive features in scattering observables when- 
ever a is increased by the factor e^^^^^22.7 [4, 5]. 

The main processes destabilizing ultracold atomic 
and molecular gases are three-body recombination, 
B-\-B-\-B^B2^B for bosons 5, and vibrational relax- 
ation, B2-\-B^B2-\-B. It is well established that the 
three-body recombination rate Ks for bosons displays 
a series of interference (a>0) and resonance (a<0) fea- 
tures as a function of a (for \a\:^ro) reflecting Efimov 
physics [2, 6, 7]. Explicitly [2, 32], 
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atomic mass. Atomic units are used here and throughout 
this Letter. While their overall behavior is determined by 
a, any comparison with a physical system also requires 
the short-range three-body parameters <l> and 77 which 
cannot, in general, be predicted from two-body physics 
alone [32]. Similar universality is also manifested by the 
a>0 vibrational relaxation rate V^ei foi" bosons [2], 
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and for FF"" ^F-^FF' collisions [20, 21] 
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where F and F' are fermions in different spin states and 
C depends on short-range physics. 

The observation of Efimov physics in ultracold gases 
has recently become a very active area [8, 9, 10, 11, 
12, 13], yet the universal behavior of three-body systems 
near a narrow Feshbach resonance remains to be fully 
understood. While specific systems have been modeled 
near a resonance [16, 17, 18, 19], no simple analytical 
expressions like Eqs. (l)-(4) for the collision rates have 
yet appeared. Instead, the vast majority of our knowl- 
edge of three-body collisions, including all of the uni- 
versal analytic expressions [2], uses only the fact that a 
Feshbach resonance changes a and ignores all other as- 
pects of the resonance. It is generally assumed that this 
approximation works well near broad resonances and less 
so near narrow ones [14] . This expectation is based on the 
low-energy expansion of the two-body 5-wave scattering 
phase shift J, 
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where sq^ 1.00624 is a universal constant and m is the 



where k is the wave number and reff is the effective range. 
Since |reff| is inversely proportional to the resonance 
width [14], the second term is negligible for broad res- 
onances, but has significant consequences for narrow res- 
onances — even at ultracold collision energies. To be a 
narrow resonance in this language requires |reff|^ro, and 
ro is on the order of 50-100 a.u. for most alkali species. A 
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typical narrow resonance for Rb with a width of 1 G cor- 
responds to reff~-200 a.u. while a quite narrow one, with 
a width of 0.1 G, has reff^-2000 a.u [14]. We should thus 
expect substantial deviations from the zero-range, ro^O, 
results in Refs. [14, 23]. 

With this Letter, we aim to expand our knowledge of 
universality near a narrow Feshbach resonance. In par- 
ticular, we study the dependence of ultracold three-body 
collision rates on reff both numerically and analytically, 
obtaining analytic expressions for and Kei similar to 
Eqs. (2) and (3), verified by essentially exact numerical 
solutions. For identical bosons, we show that in con- 
trast to [14], short-range three-body physics is important 
— even in the limit of a zero-range two-body potential. 
We also find a |reff|~^ suppression of inelastic losses for 
all three-boson processes leading to deeply bound two- 
body states. For systems involving two distinguishable 
fermions, however, we have found that large reff has the 
opposite effect: it enhances the loss rates and alters their 
universal behavior. Both the suppression and the en- 
hancement can be traced to the same modification of the 
effective three-body interaction first identified in [14]. 

Previous studies have addressed some aspects of these 
processes near a narrow Feshbach resonance and their de- 
pendence on Teff. In an especially relevant work, Petrov 
Ref. [14] calculated the recombination rate for identical 
bosons with a>0 using a zero-range potential modified 
to include rgff. By solving his model numerically, he 
found that the minima described by Eq. (1) still hold 
for a^|reff| but occur at fixed values of |reff|/a without 
any reference to a three-body parameter like <l>. Gogolin 
et al. [23] reproduced these results with a very different 
method but effectively the same physical model. Math- 
ematically, it is true that no three-body parameter is 
needed since the inclusion of reff in the zero-range model 
regularizes it [24], removing the Thomas collapse [25]. 
Physically, however, we will show that a three-body pa- 
rameter is still generally needed to represent the short- 
range three-body physics of any realistic system and that 
Teff alone is insufficient. 

To study the influence of the resonance width on the 
above collision rates, we numerically solved the three- 
body Schrodinger equation using the adiabatic hyper- 
spherical representation [7, 33]. In this representation, 
the Schrodinger equation is reduced to coupled equations 
for the hyperradial channel functions Fi,{R)^ 
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describing both bound and scattering properties of the 
system. The hyperradius R represents the overall size 
of the system, \i=ml\f?) is the three-body reduced mass, 
Wvv{R^ are the effective three-body potentials obtained 
by diagonalizing the adiabatic (flxed K) Hamiltonian, 
and Wyyi{K) are the non-adiabatic couplings governing 



transitions between channels. 

Although a Feshbach resonance is a multichannel phe- 
nomenon, we model it with a single-channel two-body in- 
teraction, since the asymptotic two-body scattering wave 
functions can be made identical for only one open chan- 
nel. In other words, the asymptotic wave function does 
not depend on how the resonance is generated [14, 15]. 
The three bodies thus interact through a pair-wise sum 
of two-body model potentials with a barrier such as 



v(t) = -L)sech^(3r/ro) + ^g-^^^^/^o-^)' 
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which has an s-wave shape resonance. In (7), the poten- 
tial depth D primarily controls a, and the barrier height 
B is adjusted to produce the desired reff. 

Figure 1 shows our numerical calculations for and 
Vre\. We generated each curve by varying a at flxed |reff |, 
corresponding to experiments tuning across a Feshbach 
resonance with a particular width. For |a|^|reff|, the 
rates retain the features predicted in Eqs. (l)-(3) as ex- 
pected. For |a|<|reff|, the rates deviate from these for- 
mulas, with for a>0 in Fig. 1(a) approaching the 
(a^ I reff I )^/^ behavior predicted in [14]. The flgure also 
shows that the rates seem to converge to a universal 
curve as |reff| increases once scaled properly. Moreover, 
the position of the flrst Eflmov feature that appears as 
|a|/|reff| increases agrees reasonably well with the zero- 
range model predictions from Refs. [14] and [23]. Note 
that the a>0 prediction is taken directly from [14], 
while we have determined the K^, a<0 prediction based 
on the location of the flrst Eflmov state in [23] and the 
Vre\ prediction based on the position of the flrst pole in 
the atom-dimer scattering length given by [14]. Recall- 
ing that reff is only about -200 a.u. for a 1 G wide 
resonance, though, we see that typical experiments will 
actually fall far from the zero-range results. To min- 
imize non-universal impacts of the two-body potential 
barrier, we have introduced a short-range hard wall in 
the Wyjj{R)^ making the short-range three-body physics 
more uniform as a function of reff. We have verifled that 
the results presented in Fig. 1 are not sensitive to the 
precise location of this hard wall. 

One of the advantages of the adiabatic hyperspher- 
ical representation is that the effective potentials of- 
fer a reasonably simple picture of the scattering pro- 
cesses [6, 7, 21]. Figure 2 shows a sketch of the effec- 
tive three-body potentials for three identical bosons with 
<^^|^efr| distilled from and verifled by numerical solutions 
of the adiabatic equation for many different combinations 
of reff, a, and v{r). The main effect of this new length 
scale, reff, is that Wyy{R) for the weakly-bound molecu- 
lar channel is modifled in the range ro<Ci^<C|reff |, taking 
the Coulomb-like form co/(2/i|reff |i^) [14] instead of the 
usual attractive R~^ Eflmov potential. Our numerical 
analysis indicates, however, that the coefficient cq is not 
universal — even changing from attractive to repulsive 
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FIG. 1: Inelastic zero-energy three-body collision rates for 
identical bosons. The vertical dashed lines indicate the pre- 
dicted position of the first Efimov feature (see [14, 23] and 
text), (a) Recombination rates for a > with one or more 
two-body s-wave bound states as indicated. The thick solid 
line is the analytical result from [14] scaled to match our data, 
(b) Recombination rates for a < 0. (c) Relaxation rates for 
a > 0. The solid lines in (b) and (c) are simply to guide 
the eye, while the dashed lines are the analytical results from 
Eqs. (9) and (8), respectively, using the a and f3 indicated. 
These were fit in the limit Ireffl^co. 



when a second two-body s-wave state is added and be- 
coming increasingly repulsive as more states are added. 
This non-universality makes it even more surprising that 
the rates in Fig. 1 are universal, which we have veri- 
fied by repeating the calculations with other model po- 
tentials [31] and with varying numbers of deeply bound 
two-body states (also shown in the figure). 

The simple picture represented by Fig. 2 can also be 
used to derive quantitative, analytic expressions for the 




FIG. 2: The schematic three-body potentials and recombina- 
tion pathways with ro^|reff|^a for a>0. 



collision rates. For the case a>0 shown in the figure, 
for instance, we calculate V^ei by considering incidence 
in the weakly bound atom-dimer channel (the lower po- 
tential in the figure) with the transition to the deeper 
two-body channels driven by non-adiabatic coupling lo- 
calized in the region R<ro. The hyperradial scattering 
wave function is found by matching the analytic solu- 
tions from each region at the boundaries a|reff| and P\a\, 
where a and p reflect the uncertainty in defining these 
boundaries and will be determined by fitting our final ex- 
pression to the numerical results. We parametrize the so- 
lution in the short-range, non-universal, region 0<R<ro 
with a complex scattering length A whose imaginary 
part accounts for inelastic transitions. In the other non- 
universal Coulomb-like region, ro<i^<a|reff|, we take the 
potential to be zero since cq can be positive or negative 
without affecting the numerical results. The relaxation 
rate is then obtained from 1/^^ —7r{l—7V)//j.k where 1Z is 
the elastic scattering probability obtained from the scat- 
tering wave function described above. The result is 
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where sinh 7^=|ImA/reff| csc(2(/:)o) sin^(<l>+(/:)o) with 
(p=so\n{l3/a)-\-(po^ tan(/?o=2so, and tan<l>=2so(a — 
ReA/|reff|)/(Q^ + ReA/|reff|). With a virtually identical 
analysis, K^^^^^ can be derived to be 



(a<o) 12\/37r^;^^ sin 2(/:?o sinh 27^ 

sin^ [so In ( I a/reff I ) + ^ + (/:'] + sinh^ rj m ' 



(9) 



Note that a and p can take on different values than for 
^ and that here tan (/9o=so/2. Similar expressions can 
be derived for other low-energy scattering observables. 

Comparison of Eqs. (8) and (9) with Eqs. (3) and (1), 
respectively, shows that the expectation that ro is re- 
placed by |reff| is indeed true but that there are also other 
modifications due to |reff|. These expressions show, for 
instance, why scaling the rates in Figs. 1(b) and 1(c) 
by |reff| produces curves with comparable magnitude: 
the factor sinh 2/^ introduces a |refr|~^ suppression. The 
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suppression in r^, which represents transitions to deeply- 
bound two-body states, is also responsible for the more 
pronounced mimina in Fig. 1(a) as |reff| increases for 
those calculations with multiple two-body bound states. 
The observation of interference minima in is thus 
more favorable near a narrow Feshbach resonance. 

Equations (8) and (9) further reveal the importance of 
the short-range three-body physics through their depen- 
dence on A in both rj and <l>. This physics is absent from 
the zero-range treatments [14, 23], so the agreement in 
Fig. 1 between our numerical results and the zero-range 
potential predictions for the position of the first Efimov 
feature is rather fortuitous. We see from the arguments 
of sin^ in our analytical expressions that A-independent 
Efimov feature positions — as predicted in [14, 23] — 
are found only in the limit ReA/lreffl^O. For the nu- 
merical examples presented in Fig. 1, Re A has a value 
comparable to ro, but this will not be true in general. 

To get a first sense of the effect of large |reff | on fermion 
collisions, we have studied the process FF'^F. For 
broad resonances |reff|— ro, the a~^"^^ suppression of the 
relaxation rate (4) originates from a repulsive barrier in 
Wyy{R) in the range ro<^R<^a [21]. When |reff|>ro, 
however, Wj^y{R) is modified in the range ro<Ci?<C|reff| 
by the emergence of a Coulomb-like potential just as 
for bosons. Therefore, the repulsive barrier in that re- 
gion is weakened, leading to enhanced vibrational re- 
laxation. Moreover, when a<c|reff|, the dependence of 
the rate on a is altered, much like the situation for 
bosons shown in Fig. 1(a). All of these effects can be 
seen in our numerical calculations shown in Fig. 3. For 
^i<keff|5 relaxation scales as (a/|reff|)~^, a much weaker 
suppression than for broad resonances. For a>|reff|, we 
find Keio^(<^/keff|)~^'^^ which recovers the expected scal- 
ing with a, but has a much larger overall magnitude 
due to the scaling with reff. This result is consistent 
with the experimental observation that ^^K2 dimers were 
found [29, 30] to have much shorter lifetimes near a nar- 
row resonance than ^Li2 dimers close to a broad reso- 
nance [27, 28]. 

To summarize, we have studied ultracold collisions of 
three identical bosons and of mixed-spin fermions near a 
narrow Feshbach resonance using the connection between 
Teff for the two-body interaction and the width of the res- 
onance. We solved the Schrodinger equation numerically 
for a wide range of parameters and showed that there is 
still universal behavior even for finite range two-body po- 
tentials. We were able to identify the key modifications 
to the three-body adiabatic hyperspherical potentials and 
thus derive analytical expressions for the rate constants. 
From these analytical expressions, we showed that short- 
range three-body physics is still important, even near a 
narrow Feshbach resonance. This result is, perhaps, un- 
fortunate for experimentalists since the positions of the 
Efimov features are, in general, still dependent on short- 
range physics, making it difficult to locate a priori a min- 
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FIG. 3: Relaxation rates for mixed-spin fermions with a>0 
and large |reff|. The relaxation rate for small |reff| is also 
plotted, showing the same scaling with a but not with |reff| 
since it is not in the universal limit. 



imum of as suggested in [14]. On the other hand, our 
analysis has shown that bosonic recombination and relax- 
ation to deeply-body two-body states is suppressed near 
a narrow resonance which might prove beneficial experi- 
mentally. Similarly, our analysis suggests that long-lived 
weakly bound FF' molecules are most easily obtained 
near a broad resonance. 

This work was supported by the National Science 
Foundation. 
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